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Abstract 

Using Jones' quadratic tangles formulas, we automate the construc- 
tion of the 4442, 3333, 3311, and 2221 spoke subfactors by finding sets 
of 1-strand jellyfish generators. The 4442 spoke subfactor is new, and 
the 3333, 3311, and 2221 spoke subfactors were previously known. 

1 Introduction 

In this paper, we construct a number of subfactors by giving explicit generators 
and relations for the corresponding planar algebras. In particular, we use 
Jones' analysis of quadratic tangles to systematically derive 'jellyfish relations', 
which are straightforward to analyze. We construct one new subfactor as well 
as three previously known subfactors. The main result of the paper is: 

Theorem 1.1. There is a '444%' subfactor at index 3 + \fh with principal 
graphs 




At this point, there are very strong number theoretic and combinatorial 
constraints on the possible principal graphs of subfactors. Indeed, when we 
find a graph that satisfies all these constraints, we reasonably suspect that 
there are actually subfactors with that principal graph. Nevertheless, the final 
step of actually constructing such subfactors remains very difficult. 

It is known that every subfactor planar algebra embeds in the graph pla- 
nar algebra of its principal graph [JP11, MW]. Thus, a standard approach 
to constructing a planar algebra is to identify some candidate elements in the 
appropriate graph planar algebra, then prove that the subalgebra they gener- 
ate is a subfactor planar algebra, with the desired principal graph. In many 
cases, it is not that difficult to obtain the candidate elements, for example, by 
solving certain polynomial equations in the graph planar algebra, or finding 
flat elements in the graph planar algebra with respect to some connection. 



The challenge is then to show that these elements generate the desired pla- 
nar algebra. The most difficult step is merely to show that they generate an 
evaluable planar algebra; that is, every closed diagram is a multiple of the 
empty diagram (equivalently, the zero box space is 1-dimensional). With this 
established, we have some subfactor planar algebra, and one wants to identify 
the resulting planar algebra. Often, as is the case here, analysis of some small 
projections in the subalgebra, and some combinatorial arguments, suffices to 
determine the principal graph. 

In [BMPS09], Bigelow-Morrison-Peters-Snyder constructed the extended 
Haagerup planar algebra, which had long been expected to exist, with principal 
graphs 



The essential insight was the jellyfish algorithm, introduced therein, which 
provides a powerful framework for proving a planar algebra is evaluable. 

Suppose we have a set of elements in a planar algebra, each a lowest weight 
rotational eigenvector, which we are thinking of as generators. A 'jellyfish 
relation' is an identity in which the left hand side is simply a single generator 
with some number of strands between it and the starred point on the boundary, 
and the right hand side is some linear combination, in every term of which every 
generator is adjacent to the starred point on the boundary. In [BMPS09], there 
were two jellyfish relations (here n = 4 corresponds to the Haagerup planar 
algebra, and n = 8 corresponds to the extended Haagerup planar algebra): 
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These are actually box jellyfish relations; it is easy enough to see that by 
expanding out the Jones- Wenzl idempotent on the left hand side in terms of 
Temperley-Lieb diagrams, and moving all non-identity terms to the right, that 
these relations become jellyfish relations as described above. See Subsection 
2.5 for more details. 

In a fc-strand jellyfish relation, the left hand side has k strands above the 
generator. Above, we have a one-strand jellyfish relation for S, the Fourier 
transform, and a two-strand jellyfish relation for S. 

A complete set of jellyfish relations is one such that by repeated application, 
we can rewrite any diagram as a linear combination of diagrams in which every 
generator is adjacent to the 'outside' starred region. We picturesquely refer to 
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this process as 'the jellyfish algorithm': we gradually float all the jellyfish to 
the surface of the ocean, possibly creating new jellyfish along the way. The pair 
of relations above is a complete set: the one-strand relation removes all the 
instances of S, then the two-strand relation allows us to float all the instances 
of S to the surface. 

It is typically easy to see that any closed diagram with all generators ad- 
jacent to the boundary is evaluable by iteratively finding an adjacent pair 
of generators which are connected by sufficiently many strands. The entire 
jellyfish algorithm is somewhat unusual amongst algorithms for simplifying a 
planar diagram; at intermediate steps, it requires making the diagram much 
more complicated. 

In this paper, we present a systematic approach to identifying jellyfish 
relations using Jones' paper analyzing quadratic tangles [Jon03]. In particular, 
we show how to find all 1-strand jellyfish relations for which the right hand 
side involves at most two generators in each term. We find that for quite a 
number of potential principal graphs, the jellyfish relations obtained in this 
way constitute a complete set. This assures us that we have constructed some 
subfactor planar algebra, and a little separate work in each case identifies the 
principal graph as the intended graph. 

It is worth noting, however, that this approach is far from uniformly suc- 
cessful! In [BMPS09], the authors needed 2-strand jellyfish relations. In- 
deed, a result of Bigelow and Penneys [BP 12] shows that having a com- 
plete set of 1-strand jellyfish relations implies that both the principal and 
dual principal graphs are spoke graphs. Thus by the triple point obstruction 
[Haa94, MPPS12], any subfactor with principal graphs beginning with a triple 
point cannot have a complete set of 1-strand jellyfish relations. We anticipate 
the results of [BP 12] giving strong constraints on subfactors whose principal 
graphs are not both spokes. 

Even with these limitations, we do have a number of interesting examples. 
We give a jellyfish presentation of the new 4442 subfactor along with the 
3311 subfactor [GdlHJ89], the 2221 subfactor (reproducing all the work of 
Han's thesis [HanlO], in an entirely automated fashion!), and one of the 3333 
subfactors previously constructed by Izumi. These graphs appear in Figure 1. 




Figure 1: The principal graphs of the four subfactors we construct in this 
paper. In each case, the dual principal graph is the same as the principal 
graph. 

Immediately after hearing about our construction of a 4442 subfactor, 
Izumi noticed that a 4442 subfactor can be constructed using a Z/3Z quo- 
tient of the 3333 subfactor. This will appear in one of his forthcoming papers. 
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Conversely, we constructed the 3333 subfactor only after hearing Izumi's al- 
ternate construction. 

In addition to giving a generators and jellyfish relations presentation of 
each of these four subfactor planar algebras, we show that the 4442, 3333, 
and 2221 subfactors are each self-dual, and moreover symmetrically self-dual. 
This essentially means that one can ignore the shading in the planar algebras, 
and thus there exist fusion categories with the same principal graphs. We will 
investigate this further in a future paper [MPP12]. 

Acknowledgements. The authors would like to thank Vaughan Jones and 
Masaki Izumi for many helpful conversations. Both authors were supported by 
DOD-DARPA grants HROOll-11-1-0001 and HR0011-12-1-0009. Scott Morri- 
son was at the Miller Institute for Basic Research and the Australian National 
University while writing this paper, and in part supported by the ARC DE- 
CRA grant DE120100232. David Penneys was also supported by NSF grant 
DMS-0856316. 

1.1 The FusionAtlas 

This paper relies on some substantial calculations. In particular, our efforts to 
find the generators in the various graph planar algebras made use of a variety 
of techniques, some ad-hoc, some approximate, and some computationally 
expensive. This paper essentially does not address that work. Instead, we 
merely present the discovered generators and verify some relatively easy facts 
about them. In particular, the proofs presented in this paper rely on the 
computer in a much weaker sense. We need to calculate certain numbers 
of the form Tr(STR), where S, T, and R are rather large martrices, and 
the computer does this for us. We also entered all the formulas derived in 
this paper into Mathematica, and had the computer automatically evaluate 
the various quantities which appear in our derivation of jellyfish relations. 
As a reader may be interested in seeing these programs, we include a brief 
instruction on finding and running these programs. 

The arXiv sources of this article contain in the code subdirectory a number 
of files, including: 

• Generators .nb, which reconstructs the generators from our terse de- 
scriptions of them in Appendix A. 

• QuadraticTangles .nb, which calculates the requisite moments of these 
generators, and performs the linear algebra necessary to derive the jel- 
lyfish relations. 

• GenerateLaTeX . nb, which typesets each subsection of Section 4 for each 
planar algebra, and many mathematical expressions in Appendices A 
and B. 
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The Mathematica notebook Generators .nb can be run by itself. The 
final cells of that notebook write the full generators to the disk; this must 
be done before running QuadraticTangles .nb. The QuadraticTangles .nb 
notebook relies on the FusionAtlas, a substantial body of code the authors 
have developed along with Emily Peters, Noah Snyder and James Tener to 
perform calculations with subfactors and fusion categories. To obtain a local 
copy, you first need to ensure that you have Mercurial, the distributed version 
control system, installed on your machine. With that, the command 

hg clone https://bitbucket.org/fusionatlas/fusionatlas 

will create a local directory called fusionatlas containing the latest version. 
In the QuadraticTangles .nb notebook, you will then need to adjust the paths 
appearing in the first input cell to ensure that your local copy is included. After 
that, running the entire notebook reproduces all the calculations described 
below. 

We invite any interested readers to contact us with questions or queries 
about the use of these notebooks or the FusionAtlas package. 

2 Background 

2.1 Graph planar algebra embedding 

A finite depth subfactor planar algebra embeds in the graph planar algebra of 
its principal graph [JP11, MW]. We begin by assuming that the desired planar 
algebra is n — 1 supertransitive, so Pfc,+ = TL k + for k — 0, . . . , n — 1 and is 
generated by the orthogonal complement P„ i+ TL n>+ of the Temperley-Lieb 
algebra in the n-box space. If this orthogonal complement has dimension k, 
we say that the planar algebra is n — 1 supertransitive of excess k. The excess 
can be read off from the principal graph; it's the sum of the squares of the 
multiplicities of the edges immediately beyond the branch point, minus one 
(or simply two less than the valence of the branch point, in the simply laced 
case) . 

To construct the subfactor planar algebras in these papers, we first find a 
connection on the desired pair of principal graphs (r,r'). We do so by first 
computing a numerical approximation to high precision, then guessing the 
exact solution using Mathematica's RootApproximant function, and finally 
verifying using exact arithmetic in a number field that the candidate exact 
solution satisfies the necessary equations. As we use this numerical method, 
we cannot say with certainty how many connections exist on these graphs. 
However, Izumi has determined all such connections, and we believe our list 
of possible connections coincides with his. As the main goal of this paper is to 
construct the new 4442 subfactor by the development of methods for deriving 
jellyfish relations, rather than proving uniqueness, we restrict our attention to 
a single connection. 
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Using this connection, it is possible to discover the rotational eigenval- 
ues of the desired generators via [IJMS11, Theorem 1.7]; in a certain gauge, 
tr([/C/*) = 2 + J^a;, where U is the 'branch matrix' of connection entries going 
through the branch point, and the sum is over the k rotational eigenvalues u, 
with multiplicity, of the lowest weight vectors in the n-box space. We then 
compute the equations for flatness and find k flat n-boxes with respect to the 
connection. (Note that these flat elements are not necessarily self-adjoint as 
in [Jon03].) 

At this point, the planar algebra generated by the flat n-boxes is some sub- 
factor planar algebra, which we denote by P„ r , and thus we have constructed 
some subfactor by [Ocn88, EK98, Jonll, MW]. However, there is important 
information about P„ r that we still need. First, we need to calculate its princi- 
pal graph (which we are hoping is the graph we started with, (r, T')). Second, 
we would like to know how to evaluate closed diagrams in the generators. 

For the spoke graphs we consider in this paper, T and T' always coincide, 
and are drawn from the set {4442, 3333, 3311, 2221} (see Figure 1), where the 
numbers refer to the length of the spokes emanating from the central vertex. 
By [BP 12], if such subfactor planar algebras were to exist with the desired 
principal graphs, they would have 1-strand jellyfish generators at depth n. We 
use Jones' quadratic tangles techniques [Jon03] to find these 1-strand jellyfish 
relations in P„ r . We then use the generators to determine that Pj" has the 
correct principal graph. 

2.2 Quadratic tangles 

In [Jon03], Jones uses quadratic tangles techniques to prove a number of for- 
mulas about planar generators in a subfactor planar algebra. In fact, many of 
the formulas there hold in more generality. 

Notation 2.1. Recall that the Fourier transform J 7 is given by 




For a rotational eigenvector S G P n ,± corresponding to an eigenvalue Us — cr|, 
we define another rotational eigenvector 5* G P njT by S = a~ s 1 J-'(S). Note that 
T{S) = a s S, so S v2 = S. 

Definition 2.2. Suppose P. is a unitary, spherical, shaded planar algebra 
with modulus 5 > 2 which is not necessarily evaluable, i.e., we don't know if 
dim(P nt ±) < oo or if dim(P 0> ±) = 1. 

A finite set 03 C P n ,+ is called a set of generators if the elements of 23 are 
self-adjoint, low-weight eigenvectors for the rotation, i.e, for all S G 53, 

• S = S\ 
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• S is uncappable, and 

• p(S) = usS for some n-th root of unity us- 

Given a set of generators 23, we get a set of dual generators 23 = {S\S G 23}. 

We say a set of generators 23 has scalar moments if Tr(R), Tt(RS), Tr(RST) 
are scalar multiples of the empty diagram in P ,+ for each R,S,T e 23. Note 
that 23 has scalar moments if and only if 23 does. 

If a set of generators 23 has scalar moments, we say it is orthonormal if for 
all S, T e 23, (S, T) = Tr (ST) = 5 S , T - 

Notation 2.3. Since 5 > 2, given a generator R e 23, the annular tangles 




UoCR) Ui(R) U 2 (R) U 3 (R) 



are a basis for !l n+ i(i2), the annular consequences of i? inside P n +i,+ [JonOl, 
Jon03]. One calculates the dual annular basis |Uj(i?) \i — 0, . . . , 2n + l} from 
(Uj(i?), Uj(-R)) = (Jjj, where the inner product is linear on the right . See 
Definition 4.2.6 of [Jon03] for an explicit formula. 

Theorem 2.4. All the formulas of §4 of [JonOS] hold in any unitary, spherical, 
shaded planar algebra with modulus S > 2 for any orthonormal set of generators 
23 with scalar moments. 

Proof. Jones explicitly restricts to an evaluable planar algebra (in fact, he says 
'subfactor,' but evaluable is the only condition we're now leaving off) before 
proving these formulas, but upon reading through the proofs, it is clear that 
the generators having scalar moments is sufficient to work in the generality we 
need here. □ 



2.3 Spherical and lopsided planar algebras 

In [MP 12], Morrison and Peters describe the spherical and lopsided conven- 
tions for planar algebras. In the spherical convention, both shaded and un- 
shaded contractible closed loops count for a multiplicative factor of 5; how- 
ever, in the lopsided convention, shaded contractible closed loops count for 1 
while unshaded contractible closed loops count for 5 2 . The main advantage 
of working in the lopsided planar algebra is that there are fewer square roots, 
so arithmetic is easier. In particular the number field in which we calculate 
is much smaller. Hence we use the lopsided planar algebra to compute the 
moments of our generators. 

The map \. p^ herical _>. prided from [ MP12 ] jg not a planar algebra map, 
but it commutes with the action of the planar operad up to certain scalars. 
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When we draw our tangles in the standard form where each input and output 
disk is a rectangle with the distinguished interval on the left and the same 
number of strings attach to the top and bottom of each rectangle, then there 
is a power of 8 for each critical point which is shaded above, and the power 
of 5 corresponds to the sign of the critical point: 



<— »> 8 



Example 2.5. We will work out the correction factor arising when commuting 
h and J 7 . Note that 



if n is even 




if n is odd. 



Hence we have 

if n is even 
1 S' 1 ^ if n is odd. 

Example 2.6. Similarly one can calculate that 

Tr(S) = ^Tr(S) = if ^ iseven 
V ; V ; \5Tr(^) if n is odd 

y ' y ' \8- Xr Yx{\S) if n is odd. 

Hence if Si, . . . , Sk G 53, we compute the moment 

8- k Tr( J'(tlS'i) • • • J^Sfc)) if n is even 



Tr(F(Si)---F(S k )) 



which implies 



Tr(5i ---Sk 



6- 1 Tr{J r {\\S 1 )---J r {\\S k )) if n is odd 



8- k Tr^Si---\\S k ) if n is even 
8~ l Tr{\Si ■ ■ ■ \\S k ) if n is odd. 



2.4 The jellyfish algorithm 

The jellyfish algorithm was invented in 
Haagerup subfactor planar algebra with 




[BMPS09] to construct the extended 
principal graphs 



t . i . i . i 
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One uses the jellyfish algorithm to evaluate closed diagrams on a set of gener- 
ators. There are two ingredients: 



'1) The generators in P n ± satisfy jellyfish relations , i.e., for each generator 
S, T, 



j(S) 




, f(S') 




can be written as linear combinations of trains , which are diagrams where 
any region meeting the distinguished interval of a generator meets the 
distinguished interval of the external disk, e.g., 



2n\ 2nT 



2nT 



r 



where Si,...,Si are generators, and T is a single Temperley-Lieb dia- 
gram. 

(Note that j(S),j(S) means the same thing as U (5'), U (5'), but we will 
use the j notation to emphasize its importance to the jellyfish algorithm.) 

(2) The generators in P Ui ±, together with the Jones-Wenzl projection f( n \ 
form an algebra under the usual multiplication 




(Note that the Mathematica package FusionAtlas multiplies in this 
order; reading from left to right in products corresponds to reading from 
bottom to top in planar composites.) 

Given these two ingredients, one can evaluate any closed diagram using the 
following two step process. 

(1) Pull all generators S to the outside of the diagram using the jellyfish 
relations, possibly getting diagrams with more S"s. 

(2) Use the algebra property to iteratively reduce the number of generators. 
Any non-zero train which is a closed diagram is either a Temperley-Lieb 
diagram or has two generators S, T connected by at least n strings, giving 
ST. 
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Section 3 is devoted to our procedure for computing the jellyfish relations 
necessary for the first part of the jellyfish algorithm, while the second part is 
rather easy 

One can see that if {A, B, f^} span a subalgebra of P n ,+, their structure 
coefficients must be given by 

„ Tr(STR) 
as ' T ~ Tr(^) ' 

and thus determined by the moments given in Appendix B. We check that the 
algebra generated by {A, B, p n >} is closed under multiplication directly in the 
graph planar algebra, in Lemma 4.1. 

Note that if we have an orthonormal set of generators, then in the no- 
tation of [Jon03], a§ T = Tr{STR) = a s R T . There are similar easy calcula- 
tions to determine the structure coefficients (3§ T of the algebra generated by 
{A, B, C P n ,-, and in the orthonormal case, (3§ T = 6^ T . 

In [BMPS09], they found 2-strand jellyfish relations , i.e., j(S) and j 2 (S) 
lie in the span of the trains of S. For spoke subfactors, we can find 1-strand 
jellyfish relations by [BP12], i.e., for each S G 53 and S G 53, j(S) lies in the 
span of the trains from 53, and j(S) lies in the span of the trains from 53. 
These 1-strand jellyfish relations are sufficient to evaluate all closed diagrams 
from our generators, and thus 53 generates some subfactor planar algebra. 

We need further arguments to prove that the resulting subfactor planar 
algebra has the desired principal graphs. It turns out that at the relatively 
low index of 3 + \/5 these arguments are easy. 



2.5 Jellyfish to box jellyfish and back again 

Given a set of jellyfish relations, we may write them in a more compact form 
in which we multiply the diagrams by a Jones- Wenzl idempotent to get rid of 
simpler diagrams. We present the arguments back and forth in the case of one 
generator, and it is clear how to generalize to the multi-generator case. 
If we know 



j(S) 




7s,s' 



n - 1 ^~ 



2n+l 



S 



Th+1 



Th + 1 



i=l 



where X G TL n+ i t+ with all strings turned down, and 



Ui(5) 



S 



, Uj(£) = @ , and U' 2n+1 (S) = Q 

i-l/ ^ \2n-i+l 2n 



(note that the U^(5')'s for i = 1, . . . , 2n + l can be obtained from the Uj-(£)'s for 
j = 1, . . . , 2n + 1 by applying suitable powers of the rotation and multiplying 



10 



by suitable powers of as), applying the Jones- Wenzl j( 2n + 2 ) to the bottom of 
the diagram gives the simpler box jellyfish relation 



f-j(S) = 



Conversely, given a box jellyfish relation of the above form, we get a jellyfish 
relation by expanding the Jones- Wenzl idempotents. 

(1) First, note that the coefficient of the identity in the Jones- Wenzl idem- 
potent j( 2n + 2 ) is 1. This term gives us a j(S) on the left hand side and 
an S o S on the right hand side. 

(2) On the left hand side, all other terms of the Jones- Wenzl either cap off 
the S, giving zero, or there is exactly one cup on the top and one cap on 
the bottom, and the cup on the top is between the 1-st and 2-nd strings 
or the (2n + l)-th and (2n + 2)-th strings. In this case, we get a scalar 
multiple of an annular consequence of S (which is not j(S)), and these 
terms can be subtracted off to the right hand side of the equation. 

(3) On the right hand side, there are a few more options. First, if any S is 
capped off, we get zero. Otherwise there is a cup between the {n + l)-th 
and (n + 2)-th strings, and we get an an S 2 which can be written as a 
linear combination of S and Temperley-Lieb diagrams. Now if there are 
more cups on top of the term in the Jones- Wenzl, the S term vanishes, 
and we get a Temperley-Lieb diagram. If there are no extra cups, then 
we are left with some Temperley-Lieb diagrams and some scalar multiple 
of an annular consequence of S (which is not j{S)). 

Hence we get a jellyfish relation. 

3 Computing jellyfish relations with quadratic 
tangles 

In Subsection 2.1, we explained how we obtained the subfactor planar algebra 
Pj" generated by the flat elements at depth n with respect to some connection 
in the graph planar algebra of V . We now describe how to calculate 1-strand 
jellyfish relations. 

The calculations in these subsections, based on the techniques from [ Jon03] , 
rely on knowing the cubic moments of the generators and the structure coef- 
ficients for the algebra generated by {A, B, f^}. This requires a computer 
calculation, but it is no more difficult than multiplying some large matrices, 




11 



with entries in a fixed number field, then taking a trace. These moments are 
given in Appendix B. 

We perform three calculations to derive the jellyfish relations. In Sub- 
section 3.1, we find those linear combinations of the quadratic tangles which 
lie in annular consequences of the generators. In Subsection 3.2, we express 
these quadratic tangles in the basis of annular consequences of our genera- 
tors. Finally in Subsection 3.3, we invert the relations found in Subsection 3.2 
to express the relevant annular consequences back in terms of the quadratic 
tangles. 

Remark 3.1. In Subsections 3.1 through 3.3, our formulas are for orthonormal 
sets of self-adjoint generators. However, the programs in the Mathematica 
notebook Quadrat i dangles .nb are slightly more general and include the 
necessary correction factors allowing us to work with orthogonal generators 
with arbitrary norms as well. 

In Subsection 3.4, we describe two ways that we verify our formulas. Since 
the computer is doing the arithmetic, we like to verify our calculations in as 
many ways that we can think of. 



3.1 Identify quadratic tangles in annular consequences 

Given our set of generators 03 = {A, B} C P n ,+, we have the dual generators 
03 = {A, B} C P n ,-, and 03 and 03 each give 4 quadratic tangles: 

{A o A, A o B, B o A, B o B} c P n +i,+ 



and 



{A o A, A o B, B o A, B o B} c P n+1 

where 

* * 



SoT = 




Since we expect these generators to give a subfactor planar algebra with 
the principal graph being the underlying graph of the graph planar algebra 
in which they were found, we hope that some linear combinations of these 
quadratic tangles lie in the space QL(A, B, 0) of annular consequences of A, B, 
and the empty diagram (the annular consequences of the empty diagram are 
the diagrams in Temperley-Lieb) , and the same for the checked generators: 

QTAC = ®(A, B, 0) n span({A o A, A o B, B o A, B o B}) 
QTAC y = Z(A, B, 0) n span({i o A, A o B , B o A, B o B}) 

(QTAC stands for "quadratic tangles in annular consequences"). 
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Remark 3.2. In our notation 2t(A, B, 0) includes Temper ley-Lieb. When we 
refer to the annular consequences of A, B only, we will call this space %{A, B) 
(which is what is referred to simply by 9 in [Jon03]). 

Example 3.3. Starting with our 2 generators found in the graph planar algebras 
of the 4442 or 3333 principal graphs, since 4442 and 3333 each have annular 
multiplicites *22, we hope that dim(QTAC) = dim(QX4C v ) = 2, i.e., the 
quadratic tangles are as linearly independent as possible. 

Example 3.4. Starting with our 2 generators found in the graph planar algebra 
of the 3311 principal graph, since 3311 has annular multiplicites *20, we hope 
that dim(QTAC) = dim(QX4C v ) =4. 

Example 3.5. Starting with our 2 generators found in the graph planar algebra 
of the 2221 principal graph, since 2221 has annular multiplicites *21, we hope 
that dim(QTAC) = dim(QX4C v ) = 3. 

We use formulas from [Jon03] to calculate bases for QTAC and QTAC y . 
We first describe how to find a basis for QTAC. 

First, we calculate the 4x4 matrix of inner products modulo the annular 
consequences of A, B: 

(([SoT-P% (AB) (SoT)],PoQ}) 

The inner products are given by the following formulas, where the second 
comes from Proposition 4.4.2 of [Jon03]: 

(S o T, P o Q) = Tr(QT) Tr(SP) 

M 

(P%(A,B)S o T, P o Q) = J2 ttU (^¥4 Q + VT<T^v P b^b P R Q ) (c^ 1 + [2n + 2]) 

+ (-l) n+ V* (*Q*paWb p R Q + a T W^bf4 Q ) (^r> + 1]) }, 

where a s R T = Tr{RST), b s R T = Tr (RST), and W R = q 2n+2 + q~ 2n+2 -uj r -uj r 1 . 
Here, q > 1 such that 5 = [2] = q + q' 1 . 

Taking a basis for the nullspace of this matrix gives us a basis for QTAC. 

Remark 3.6. To calculate a basis for QTAC y , one passes to the dual planar 
algebra of the graph planar algebra (the same planar algebra, but with the 
shading reversed), and uses the same formulas above. This amounts to switch- 
ing a s R and b s R ' , since S and S have the same chirality for each S G 55. 

3.2 Find the jellyfish matrices 

Now we want to write each basis element of QTAC or QTAC y as a linear 
combination of the L)i(R), Uj(-R)'s for R G 53, R G 03 respectively. We describe 
the process for the basis elements of QTAC, and the checked versions are again 



13 



computed by passing to the dual as in Remark 3.6. Using Proposition 4.4. l.i 
in [Jon03], we have 



P&(A,B)(S oT) = J2 ° n R a s R T G n+1 (R) + a^a s b s R T (R). 



(3.1) 



Re's 



We can express U n+ i(i?), Uo(-R) (in the dual annular basis) in terms of the 
annular basis Ui(R) for R G 03, using the formulas from Proposition 4.2.9 of 
[Jon03]: 



COR) 
n+1 {R) 



1 



W R 
1 



([2n + 2] U (R) + ((-a) n+1 + {-a)- n - l )[n + 1] U n+1 (R) + X) 
([2n + 2] U n+1 (i?) + ((-a) n+1 + (-a)-™- 1 )^ + 1] U (R) + Y) 



where X,Y e %(R) = span {Ui(i2)|i ^ 0,n + l}. 

Remark 3.7. Our calculations, available bundled with the arXiv sources for 
this paper in the Mathematica notebook QuadraticTangles .nb, are slightly 
more complicated; we don't actually use the formulas from Proposition 4.2.9 
above, but instead directly compute the change of basis matrix from the dual 
annular basis to the annular basis as an extra check of these formulas. The 
change of basis matrix is computed as follows. 

First, letting U and U be the column vectors corresponding to the basis 
elements Uj(i?) and Ui(R), there is some matrix V G M 2n+ 2(C) such that 
VU = U. From the formulas (Uj(i2), Uj(R)) = 5 i; j and 



W, 



5 
a 




±i 

R 



if % = j 
if j = i ± 1 
else, 



we have V = W~ l (remember that the inner product is linear on the right). 

Putting it all together, we can express v in our basis of QTAC as a lin- 
ear combination of U (A), U (B), \J n+ i(A), U n+ %(B), plus another element in 
a(A,£,0) = TL n+1) ± © span({U i ( J R)| J R G 93 and i ^ 0, n + 1}). 



E 

S,Te<8 



7s,t' 



5RT 




Tn + 1 



Tn + 1 




where Z G 5, 0) and 



S 1 if n is even 
S if n is odd 



S 



Vn 



S if n is even 
S if n is odd. 
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Given such an equation, we can multiply by a Jones- Wenzl idempotent in 
two ways to find the following relations, from which we will derive the desired 
box jellyfish relations: 

(1) We can multiply by /( 2 ™+ 2 ) on the bottom to isolate the U (A),U (B): 



E 



1st 



Tn + l Tn+1 



y(2n+2) 



= Ha 



A 



\2n 



j(2n+2) 



B 



\2n 



j(2n+2) 



"i 1 r 



i 1 — n 1 r 

as any cap on top of /( 2n + 2 ) gives zero. 
(2) We can bend /( 2n + 2 ) around the top to isolate the \J n+1 (A), U n+1 (B): 





+ v B 




which is equivalent to 



\q . >,\" — '/•/• n) 



E l-n n-l j" + 1 
1S,T(T S a T I ' w „. 




n+ 1 



j(2n+2) 



V A 



y(2n+2) 

i 1 r 



+ VB 



)Vjj 



2n 



Remark 3.8. The second relation above is actually superfluous. It suffices 
to consider the relations above of the first type for both 03 and 23. On the 
other hand, the computer is doing the arithmetic, so we prefer to get a nice 
consistency check on all our formulas with little extra work. 

Notation 3.9. For S,T G 93,23, we use the notation 



s 



\2n 



f-j(s) 



s 



\2n 



j(2n+2) 



, and / • (S o T) 



sY^Ht 



ra+l n+1 



j{2n+2) 



Definition 3.10. Suppose {v± : . . . , Vk} C QTAC, {vi, . . . ,i)k} C QTAC y are 
bases (in our examples, 2 < k < 4). First, we use the above method to 
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calculate coefficients fj, A , /i B , fi A) fi l B so that 

«i = E ^ 5 ° T ) = ^[f ■ + /£[/ ■ 3(B)] and 

S,Te<B 

*« = E 75,T/(5°T)=^[/-j(^)]+/i i B[/-j(5)]. 



s,Te<8 



The jellyfish matrices are the matrices J, J whose z-th rows are (^i A /x^) , 
As) respectively. 

The quadratic tangles matrices are the matrices K, K whose 2-th rows are 

{Ya,a Ja,b,Jb,a>7b,b), {Ya,a 1a,b,1b,a^b,b) respectively. 
Note that 

ff(AoA)\ 



K 



f(AoB) 
f(BoA) 
\f(BoB)J 

and similarly for the checked version 



J 



f-j(B) 



3.3 Invert the jellyfish matrices to get box jellyfish re- 
lations 

Given the jellyfish matrices J, J, we check if they have rank 2. If they do (and 
we know they should by [BP12]), then we find a left inverse by the formula 

J- 1 = (J* J)- 1 J* 

(and similarly for J), since J* J G M 2 (C) has rank 2. We then use J -1 , J -1 to 
get the box jellyfish relations 



f-j(B) 



J~ X K 



ff(AoA)\ 
f(AoB) 
f(B o A) 

\f(BoB)J 



which express / • j(A),f ■ j(B) as linear combinations of quadratic tangles, 
and similarly for the checked versions. 



3.4 Checking our formulas 

As emphasized in Remarks 3.7 and 3.8, since the computer is doing the arith- 
metic, we like to check our formulas in as many ways as we can think of. We 
perform two extra checks of the formulas obtained from the above calculations. 

(1) Given our set of generators 03 in the graph planar algebra, we can com- 
pute the annular bases for 3(A), Q(B) and the quadratic tangles S o T 
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for S, T G 03 directly from the graph planar algebra. We can then use 
numerical linear algebra to compute an approximate basis for QTAC, 
and we can compare these results with those obtained in Subsection 3.1. 
We do a similar check for 03 and QTAC V . 

Similarly, we can compute numerical approximations for the coefficients 
of the quadratic tangles in annular consequences with respect to the 
annular basis directly in the graph planar algebra. We compare these 
numbers with those computed from Equation (3.1) using Remark 3.7. 

Finally, we can numerically find the jellyfish formulas J~ l K, J~ l K di- 
rectly from the generators in the graph planar algebra and compare them 
with the J~ 1 K, J^K computed in Subsection 3.2. 

All of these checks are carried out in Quadrat i dangles .nb, in the "QT 
Direct" sections for the graphs 3333, 3311, and 2221. 

(2) It would be beneficial to check the actual jellyfish formulas directly in the 
graph planar algebra. Since the Jones- Wenzl idempotent, written in the 
graph planar algebra, is expensive to compute, this is only feasible for 
the smallest graph 2221, and even then, we need a clever trick introduced 
in [BMPS09] by Stephen Bigelow. First, note that 



f(n+l) 



n + 1 



71+1 



2n+2) 



which implies that for S, T £ 03 and all 7, 



f(SoT) 



sY^Ht 



I n + 1 I n + 1 I n + 1 | n + 1 



j(2n+2) 



T 



T 



j(2n+2) 



T 



T 



7 



./ 



(n+1) 



n + 1 n+1 



j{2n+2) 



Note that ST — T^jjpr /^ is uncappable as it lies in span{A, B}. Hence 
applying any two caps which do not enclose * to S o T — r^^j 
gives zero, and therefore f(S o T) is equal to 




n + 1 Jn + 1 * 



Tn + l \n + l 



Tr (ST) 



(2n+2) 
cup 



[n+l][n + 2] 



f(n+l) 



n+1 n+1 



fx 



(2n+2) 
cup 



where is the sum of all terms in the Jones- Wenzl j( 2 ™+ 2 ) with 

exactly one cup on the top and one cup on the bottom. A formula 
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for the 1-cup Jones- Wenzl in terms of Temperley-Lieb diagrams can be 
deduced easily from [Rez07, Mor]: 



J l—cup 



fc-2 

E 

a=0 



+ 



E 



a + 1] [k - a - 1] 

[fc] 









a 











W/ / 












(■ 


a 




b 


, + y 


a 




6 



a+fe+c=fc-2 
c>0 

Note that the number of terms in the 1-cup Jones- Wenzl grows linearly, 
whereas the number of terms in the Jones- Wenzl idempotent, given by 
Catalan numbers, grows exponentially 

Using this trick, we numerically check the jellyfish formulas calculated 
in Section 4.4 in the graph planar algebra of 2221. 



4 Generators and relations 

We now have a subsection for each of our subfactor planar algebras. The three 
lemmas in each section show the results of the calculations described above. 
The proofs are simply substituting the appropriate quantities (moments, chi- 
ralities, etc.) into the formulas above. You can verify all these calculations 
using the Mathematica notebooks included with the arXiv sources of this 
paper. 

Throughout, the notation Xa„,...,a denotes the root of the polynomial Yli a i x% 
which is closest to the approximate real number z. (The digits of precision of 
z are in each case chosen so that this unambiguously identifies the root.) Thus 
for example A^'H o 8 -864 o 81 denotes the root of 1024a; 4 — 864x 2 + 81 which is 
closest to 0.3278. ' 

Lemma 4.1. For each of the graphs T = 4442,3333,3311 or 2221, the ele- 
ments A, B and in the (n,+)-box space of the graph planar algebra are 
closed under multiplication, and their structure coefficients 

ST = aj jT A + a B StT B + a{ T f {n) 

are given the following ratio of moments: 

R _ Tr(STR) 
as ' T ~ Tr(i?2) • 

A similar result holds for the elements A, B and f^ in the (n, —)-box space. 

Proof. The program Verif yClosedUnderMultiplication in the Mathematica 
notebook QuadraticTangles .nb verifies that the algebra generated by the set 
{f^ n \ A, B} is closed under multiplication directly in the graph planar algebra. 
Once we know this, the formula claimed for the structure coefficients follows 
by taking inner products. □ 
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4.1 4442 

Lemma 4.2. The linear combinations 



K 



(AoA\ 
AoB 
Bo A 

\BoBj 



and K 



(A 
A 
B 



oA\ 

oB 

oA 

oBj 



lie in annular consequences, where 



K 




x (0.809+0.588i) 
A l,-1, 1,-1,1 



and 



K = 




x (0.809+0.588*) 
A l -1,1,-1,1 



-r 
o 



Lemma 4.3. In particular, we have 



K 



ff(AoA)\ 
f(AoB) 
f(B o A) 

\f(BoB)J 



J 



■3(A) 
J{B) 



and K 



(f(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



where 



J 



A (7.275) 

' A 109,0,-5770,0,25 



A 



(6.745+2.19H) 
11881, 
0, 

-966285, 
0, 

30007665, 
0, 

1366875, 
0, 

164025 



x (-7.275) \ 
A 109,0,-5770,0,25 ' 

x (6.745+2.191i) 
A 11881, 



0. 



-966285, 
0, 

30007665, 
0, 

1366875, 
0, 

164025 



and 



J 



A (7.275) 

' ^109,0,-5770,0,25 

x (6.745+2.191i) 
A 11881, 



-966285, 
0, 

30007665, 
0, 

1366875, 
0, 

164025 



A 



(-7.275) > 
109,0,-5770,0,25 

x (6.745+2.191i) 
A 11881, 
0, 

-966285, 
0, 

30007665, 
0, 

1366875, 
0, 

164025 
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Lemma 4.4. The elements A and B satisfy the box jellyfish relations 



J -3(B) 



= J~ 1 K 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



and 



j{B) 



ff(AoA)\ 
f(AoB) 
f(B o A) 

\f(B o B)J 



where 



J- X K 



/, (0.06872) 
/ A 400,0, -23080,0,109 



A 



, (-0.06872) 
A 400,0, -23080,0, 109 



(0.067054-0.021787i) 
41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 



, (0.067054-0.021787i) 
A 41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 



, (0.067054+0.021787i) 
A 41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 

, (0.067054+0.021787i) 
A 41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 



,(-0.06872) \ 
A 400,0, -23080,0,109 » 



,(0.06872) 

A 400,0, -23080,0, 109 



and 



A (0.06872) 
400.0.-2.' 



400,0,-23080,0,109 



A 



(-0.06872) 
400,0,-23080,0,109 



A 



(0.067054-0.021787i) 



A 



41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 

(0.067054-0.021787*) 

41990400, 

0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 



, (0.067054+0.021787i) 
^41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 

, (0.067054+0. 021787i) 
A 41990400, 
0, 

87480000, 
0, 

480122640, 
0, 

-3865140, 
0, 

11881 



A 



(-0.06872) \ 
400,0,-23080,0,109 » 



A 



(0.06872) 

400,0,-23080,0,109 



4.2 3333 

Lemma 4.5. The linear combinations 

/AoA\ 

AoB 
Bo A 
\BoBJ 



K 



(AoA\ 
AoB 
Bo A 

\BoBl 



and K 
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lie in annular consequences, where 



K 



1 i (-3 - V5) 
Oil 



and 



K 



1 i (-3 - y/5) 
Oil 



Lemma 4.6. In particular, we have 



K 



(f{A 
f(A- 
f(B 

\f(B 



A)\ 

B) 
■A) 



1 I and K 

\f-j(B)J 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



r{f-M)\ 

\f-mj' 



where 



J 



§ (-3 - 3x/5) 



(-2 - V5) 
|(1-V5) 



and 



J = 



A 



(-1.102) 

1024,0,-1344,0,121 



A 



(-0.3278) 
1024,0,-864,0,81 



A 



x (-0.2860) 

^1024,0,-96,0,1 

(1.102) 

1024,0,-1344,0,121 



Lemma 4.7. TTie elements A and B satisfy the box jellyfish relations 



f-j(B)J 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



and 



f-j(B)) 



ff(A 
f(A- 
f(B 

\f(B 



A) \ 
B) 
•A) 

B) J 



where 



(1-VS) 

1(3-3^5) |(2-V5) |(2-V5) J (l + VS) 



5-2) K-1-V5) 
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and 



(-0.8422) 
X 64,0, -216,0,121 

(-0.2504) 
X 64,0, -1296,0,81 



x (-0.2185) 
A 64,0,-24,0,l 

x (0.8422) 

A 64,0,-216,0,121 



x (-0.2185) 
A 64,0, -24,0,1 

, (0.8422) 

^64,0,-216,0,121 



A 



(0.7349) 
5184,0,-3024,0,121 
(0.2185) 



A 



64,0,-24,0,1 
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4.3 3311 

Lemma 4.8. The linear combinations 



K 



(AoA\ 
AoB 
Bo A 
\BoBj 

lie in annular consequences, where 



K 



and K 



and 



K 



(AoA\ 
AoB 
Bo A 

\BoBj 



n 








o\ 





1 














1 





\o 








1/ 


/I 








o\ 





1 














1 





\o 








1/ 



Lemma 4.9. In particular, we have 



K 



where 



ff(AoA)\ 
f(AoB) 
f(B o A) 

\f(BoB)J 



and 



J = 



J 



J 



f-j(B) 





(-0.7113i) 
144,0,312,0,121 

x (0.7113i) 
A 144,0,312,0,121 



A 



and K 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



i(-7-3V3)\ 
l(l-V3) 
l(l-V3) 

^(7^3-9) y 



J 



(-0.1888) 
7776,0,-3672,0,121 



A 



x (-0.2832i) 
A 1536,0,1632,0,121 

x (0.2832i) 
^1536,0,1632,0,121 

(0.860) 

24576,0,-2414976,0,1771561 



A 



(-0.0621864) \ 
58806,0,-486,0,1 

x (0.4953i) 
A 864,0,216,0,l 

x (-0.4953i) 
A 864,0,216,0,l 

(0.2832) 

1,-1632,0,121/ 



f-m 

f-j(B) 
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Lemma 4.10. The elements A and B satisfy the box jellyfish relations 



f-j(A) 
f-j(B) 



= J~ 1 K 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



and 



f-j(B) 



= J~ 1 K 



ff(AoA)\ 
f(AoB) 
f(B o A) 

\f(BoB)J 



where 



J- 1 K = 



-137208-77672^ 
172163 



, (0.7029i) 
A 121,0,78,0,9 

-33708-32332^ 
172163 



x (-0.7029i) 
^121,0,78,0,9 

-33708-32332^ 
172163 



52626+80142^ / 
172163 / 



and 



J- l K 



A (-0.205052) 
I A 1559184260929, 
0, 

-65569974336, 
0, 

497664 



A 



(-0.11725) 

1559184260929, 

0, 

-23898282768, 
0, 

33872256 



x (0.27984i) 
A 3872,0,648,0,27 



x (-0.850i) 
A 512,0,4896,0,3267 



x (-0.27984i) 
A 3872,0,648,0,27 



A 



(0.850i) 

512,0,4896,0,3267 



-viu.yootoj \ 
^1559184260929, 



(0.93378) 
0, 

-1362157739172, 
0, 

2305430424 

, (0.53393) 
^3118368521858, 
0, 

-1989820501362, 
0, 

313826716467 



/ 



4.4 2221 



Lemma 4.11. The linear combinations 



K 



(AoA\ 
AoB 
Bo A 
\BoBj 

lie in annular consequences, where 

(\ 
1 



K 



1 



and K 



(AoA\ 
AoB 
Bo A 

\BoBj 



^ (-23- 7V21) \ 

x (1. 050-1. 818i) 
A 625, 
0, 

2300, 
0, 

10464, 
0, 

-7360, 
0, 

6400 

x (1. 050+1. 818i) 
A 625, 
0, 

2300, 
0, 

10464, 
0, 

-7360, 
0, 

6400 
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and 



(\ 
1 



K 



1 



i(-23-7V5l)\ 

> (1. 050-1. 818i) 
A 625, 
0, 

2300, 
0, 

10464, 
0, 

-7360, 
0, 

6400 



A 



(1. 050+1. 818i) 
625, 
0, 

2300, 
0, 

10464, 
0, 

-7360, 
0, 

6400 



Lemma 4.12. In particular, we have 



K 



f(AoB) 
f(B o A) 
\f(BoB)J 



J 



f-j(B) 



and K 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



where 



J 



(\ (-6 - V21) A 



A 



A 



(1.680-4.996i) 

2025, 

0, 

90855, 
0, 

1616571, 
0, 

931200, 
0, 

160000 

(1.680+4.996i) 

2025, 

0, 

90855, 
0, 

1616571, 
0, 

931200, 
0, 

160000 



(1.326) 

225,0,-393,0,-5 



\ 



x (-0.6319-1.0945i) 
A 9,9,12,-3,l 



A 



(-0.6319+1.0945i) 
9,9,12,-3,1 
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and 



J 



(\ (-6 - V2l) A^2 5 3 o!-393,0 -5 \ 

x (1.680-4.996*) x (-0.6319-1.0945*) 
A 2025, A 9,9,12,-3,l 
0, 

90855, 
0, 

1616571, 
0, 

931200, 
0, 

160000 

x (1.680+4.996*) x (-0.6319+1.0945?) 
A 2025, A 9,9,12,-3,l 
0, 

90855, 
0, 

1616571, 
0, 

931200, 
0, 

160000 



Lemma 4.13. The elements A and B satisfy the box jellyfish relations 



f-j(B) 



J- l K 



ff(AoA)\ 
f(AoB) 
f(BoA) 

\f(BoB)J 



and 



j{B) 



ff(AoA)\ 
f(AoB) 
f(B o A) 

\f(BoB)J 



where 



(k (33 



8V2T) A 



A 



(0.3277) 

2601,0,885,0,-125 



A 



(0.034194+0.063236*) 
4228250625, 
0, 

18810887175, 
0, 

36065983311, 
0, 

203997780, 
0, 

960400 



(-0.1561+0.1682*) 



2601,4896,18981,5700,925 /v 2601,4896,18981, 5700,925 



A 



x (0.034194-0. 063236*) 
^4228250625, 
0, 

18810887175, 
0, 

36065983311, 
0, 

203997780, 
0, 

960400 
(-0.1561-0. 1682*) 



659V21-2049 
2550 



x (-0.07717) 
A 65025,0,1149168,0, 



-6845 



and 



(33 - 8x721) 



A 



(0.3277) 

2601,0,885,0,-125 



A 



A 



A 



(0.034194+0.063236*) 
4228250625, 
0, 

18810887175, 
0, 

36065983311, 
0, 

203997780, 
0, 

960400 
(-0.1561+0.1682i) 

2601,4896,18981,5700,925 /v 2601,4896,18981, 5700,925 



A. 



(0.034194-0.063236*) 
4228250625, 
0, 

18810887175, 
0, 

36065983311, 
0, 

203997780, 
0, 

960400 

-0.1561-0.1682*) 



659V21-2049 
2550 



x (-0.07717) 
A 65025,0,1149168,0, 



-6845 
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5 Self-duality and calculating principal graphs 



We now know that our elements A,B£ VA(T) generate an evaluable planar 
subalgebra P, r , and hence a subf actor planar algebra. We now show that the 
principal graphs of the Pj" planar algebra are, in fact, (T, T). 

5.1 Self duality 

In this subsection, we show that P 4442 , P 3333 and P 2221 are self-dual, i.e., there 
is a planar algebra isomorphism P, — > P^ where P. v is the dual planar 
algebra obtained from P. by reversing the shading. Note that this means for 
all k, there is a map &k,±'- Pk,± Pk± — Pk,^-, an d these maps commute with 
the action of the planar operad. 

In fact, these three subfactor planar algebras are more than self-dual; they 
are symmetrically self-dual , i.e., for every n, $ niT o $ n ± = l n> ±. Hence by 
[MPP12], we can lift the shading to get fantastic planar algebras, i.e., un- 
shaded, spherical, evaluable C*-planar algebras. 

Given a fantastic planar algebra P,, we have an associated rigid C*-tensor 
category Cp. whose objects are the projections of P, and a morphism in 
Homc(p — y q) is an element x G P, such that x = pxq (see [MPS 10] for more 
details). Note further that Cp. is generated by a single self-dual object X (the 
strand), and Cp. is Z/2-graded. The fusion graph with respect to X is exactly 
the principal graph of P.. If the fusion graph is finite, then Cp. is a unitary 
fusion category. 

Hence the subfactor planar algebras P 4442 , p 3333 ; p222i gj ve r j ge ^ 2^2- 
graded unitary fusion categories generated by a single self-dual object with 
fusion graphs 4442,3333,2221 respectively. Note that a fusion category with 
fusion graph 2221 has previously been constructed by Ostrik in the appendix 
to [CMS11]. 

Theorem 5.1. The map P 44 ^ 2 P 44 f 2 exchanging A •<->■ A and B -H- B gives 
a symmetric self-duality of planar algebras. 

Proof. By the symmetry of the moments in Appendix B.l, the map clearly 
preserves the moments. Hence the box jellyfish relations of Subsection 4.1 are 
preserved under $. Moreover, the structure coefficients in the algebra P^^± are 
also preserved, so by the jellyfish algorithm the map preserves the evaluation 
of all closed diagrams so is an isomorphism. □ 

Theorem 5.2. The map $ 3j+ : P|^ 3 -> P 3333 by 



U) " M \b) whm M = k 

\ 4 2 
gives a symmetric self-duality of planar algebras. 
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Remark 5.3. It might be possible to choose generators so that (A, B) h-> (A, B) 
is already a symmetric self-duality of planar algebras, but it seems that one 
would have to work in a larger number field for this to be possible. 

Proof. One can easily verify that this map preserves the moments given in 
Appendix B.2, so the argument from the proof of Theorem 5. 1 applies. (In fact, 
we do this verification in the Mathematica notebook Quadrat i dangles .nb.) 
Finally, note that M 2 = 1, so $ n;T o $ n ± = l„ j± . □ 

Theorem 5.4. The map P 2221 -H- P 2221 swapping A <-> A and B <-> B gives a 
symmetric self- duality of planar algebras. 

Proof. Similar to the proof of Theorem 5.1. □ 

Note, however, that the result is not true for 3311: 
Theorem 5.5. P 33U is not self-dual. 

Proof. Suppose there were a self-duality $. Since $ is a map of planar algebras, 
$ must preserve the Temperley-Lieb planar subalgebras TLk t ± and also the 
low weight spaces for the rotation in Pl±. Since the rotational eigenvalues are 
uja = —1 and u>b = 1, we know that $(A) = fiA for some /i G C x . However, 
this clearly violates 

Tr(A 3 ) = J- (-6 - 4^) and Tr(i 3 ) = 

(see Appendix B.3). □ 

5.2 Identifying principal graphs 
Theorem 5.6. The principal graphs of P^ 2 are 




Proof. First, note that the modulus of the given subfactor planar algebra is 
y/3 + \/5 ~ 2.28825. Let T denote the principal graph. Using the moments 
calculated from our generators in Appendix B.l, we find all solutions to the 
equation 

(aA + bB + cf {5) ) 2 = aA + bB + c/ (5) 
to get that the minimal projections one past the branch point are given by the 
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equations: 



x (0.055601) 

A 332150625,0, -26222676750,0,254210027400,0,-847113120,0,190096 
> (-0.015556) 

A 332150625,0,-26222676750,0,254210027400,0, -847113120,0,190096 



A 



(a, b, c) = < 



(-0.015556) 

332150625,0,-26222676750,0,254210027400,0,-847113120,0,190096 
^ (0.055601) 

A 332150625,0,-26222676750,0,254210027400,0, -847113120,0,190096 



J(-5 + 3V5)) 



' (-0.040045) 
v A 2025,0,-271890,0,436 

x (-0.040045) 
A 2025,0, -271890,0,436 

g(i5-6V5)). 

Note that the dimensions of these projections agree with the dimensions of 
the vertices on the 4442 graph, so all the arms on the principal graph T must 
continue. Two branches cannot merge since 

=s 2.33743, 
=s 2.31725, 



no branch can split since 
and no branch can have a double edge since 



2.4761, 



and all of these numbers are already too large. Hence all branches continue 
simply. By counting Frobenius- Perron dimensions, one arm must stop, but 
the other two must continue. Again, the branches cannot merge or split as 

2.32033, 
2.29079, and 
2.41976. 

Again, both branches must continue simply. Once more by counting dimen- 
sions, the two remaining arms must continue, and again they cannot merge or 
split as 

2.30231, 
2.29193, and 
2.37309. 
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We conclude by counting dimensions again that T is 4442. 

The dual principal graph is also 4442 by Theorem 5.1. 

To determine the dual data, we run the program FindGraphPartners on 
the graph 




which determines all possible pairs of principal graphs and dual data for which 
one of the graphs is 4442. The only possibilities for which the principal and 
dual principal graph are both 4442 are 



and 



We rule out the second pair by noting that the dimension 1 bimodules form a 
group, and thus the vertices at the far right must be dual to each other. □ 

If we already knew P 4442 was finite depth, we could give an alternative 
argument identifying the principal graph as 4442 as follows. We omit details, 
as this is redundant with the argument in Theorem 5.6. 

Theorem 5.7. The only finite depth subf actor principal graphs with index 
3 + \fh starting like J^lll are 444%- 

Proof. If any of the edges above the quadruple point end immediately, that 
vertex has dimension which is not an algebraic integer (a root of 81 — 126x 2 + 
4a; 4 ). Otherwise, graph enumeration in the style of [MS10] shows that two of 
the legs end, with lengths 2, 2 or with lengths 2, 3. In second case, there's 
a dimension which isn't an algebraic integer. In the other case, the graph 
contains 4422. (There are two possibilities for the dual data.) We can look 
for connections on this graph (without assuming that it ends), and show that 
there are no bi-unitary connections. □ 

Theorem 5.8. The principal graphs o/P. 3333 are 



Proof. The proof is similar to that of Theorem 5.6. Again, the modulus is 
\/3 + V5 ~ 2.28825, and we find that the minimal projections one past the 
branch are given by a A + bB + cf^ where 



(a,b,c) 




,6( 3 -^)4( 3+ ^)-£)' 
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Since Tr(/^) = 6 + 3\/5, all the minimal projections have trace 2 + a/5, which 
agree with the Frobenius-Perron dimensions of the vertices of 



at depth 4. Once again, counting Frobenius-Perron dimensions and noting 
that 




yields the result. 

The dual principal graph is also 3333 by Theorem 5.2. 

To determine the dual data, we run FindGraphPartners as in the proof of 
Theorem 5.6 on the graph 



The only possibilities for which the principal and dual principal graph are both 
3333 are 



and 



Since uja = = 1 (see Appendix A. 2), the two click rotation p must be the 
identity on P| 3 ^ 3 0TL 4j+ , and thus all vertices at depth 4 must be self-dual. □ 

Theorem 5.9. The principal graphs o/P 3311 are 



Proof. Similar to the proofs of Theorems 5.6 and 5.8. The modulus is a/3 + \/3 
2.17533, and the minimal projections are given by aA + bB + cf^ where 



(a, b, c) 



(-^(4-3^)^(1 + 2^ 

(|(3-V3),i(7 + 3V3),i(5-V3)) 

(1 (-3 + 3a/3) , i (-15 + 3a/3) , n (5 - a/3)) 
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Since Tr(/^) = 4 + 3\/3, the traces of the minimal projections are 2 + a/3, 1 + 
\/3, 1 + \/3 respectively, which agree with the Frobenius-Perron dimensions of 
the vertices of 



at depth 4 reading from bottom to top. Finally, 



a 2.33441, 
a 2.23607 and 
a 2.47485. 

Hence the principal graph is 3311. 

For the dual graph, note that the minimal projections in P\_ are given by 
a A + bB + c/^ where 

(o, ^2 (27 -V3), i (1 + 2V3)^j 
l-lyfl (3-V3), -^(27- V3), ^t(5- V^)) 
^ (l^l (3-V3), (27- v^), i (5 - v^)) 



(a,6,c) 



Since Tr(/( 4 )) = Tr(/( 4 )), the same argument as above applies, and the dual 
principal graph is 3311. 

To determine the dual data, we note there are only two possibilities for 
each graph: either the singly valent vertices at depth 4 are self-dual or they 
are dual to each other. Since uja = — 1 and cub — 1 (see Subsection A. 3), we 
know that p 2 = 1 on P|+ n TL^ + . Hence the singly valent vertices must be 
self dual. □ 

Theorem 5.10. The principal graphs of P 2221 are 



Proof. Similar to the proof of Theorems 5.6, 5.8, and 5.9. The modulus is 
y(5 + y / 2T)/2 ~ 2.1889, and the minimal projections are given by aA + bB + 
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c/( 3 ) where 



(a, b, c) — < 
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483 - 103V21 - 5W6 -5283 + 1153V21 



33 - 7V21 - W6 -87 + 19V21 



483 - 103^ + 5-y/6 (-5283 + 1153v^l) 



{ \ 



3 (23 + 5\/2T 



103 i 
v^T) ' \ 150 50 ' 5 V J 



\ 



J 



Since Tr(/( 3 )) = v19 + 4a/21, the traces of the minimal projections are 

| + \Z|"+ ^F> v ^ respectively, which agree with the Frobenius-Perron 
dimensions of the vertices of 



at depth 3 reading from bottom to top. Finally, 




2.32437, 



2.22158, and 



2.46991. 



The dual principal graph is also 2221 by Theorem 5.4. 

Finally, since the dimension one bimodules form a group, the dual data 
must be as claimed. □ 



A Generators in the graph planar algebra 

Specifying an element of the graph planar algebra of a large graph can be a 
somewhat cumbersome process; the element is a function on loops of a certain 



32 



length on the graph, and we need to specify each value. Fortunately, if the 
element is a lowest weight vector, this can be significantly abbreviated. 

Throughout this section, we assume that T is a spoke graph with m + 1 
arms, and the initial arm is at least as long as any of the other arms. (This 
obviously holds for the graphs in which we are interested.) 

Lemma A.l. 

(1) A lowest weight vector in A e G(T) n>+ is determined by its values on 
'collapsed' loops which stay within distance one of the central vertex. 

(2) In fact, it is determined by its values on such loops which never enter 
one of the spokes of our choice. 

(3) Further, the value on any loop which has more than 2k + 1 consecutive 
vertices which either lie on a particular arm of the graph of length k — 1 
or are the central vertex is zero. 

Proof. We'll work in the spherical graph planar algebra, as it is somewhat 
easier to state the requisite formulas there. Obviously the lemma holds in the 
spherical planar algebra if and only if it holds in the lopsided graph planar 
algebra. 

Call the central vertex of T c. We'll write ||7|| = £V d{ r y{i), c). For a 
collapsed loop, ||7|| = n/2, while for any other loop ||7|| > n/2. We'll show 
that for any non-collapsed loop 7, if A is a lowest weight vector, then ^(7) is 
determined by the value of A on loops of strictly smaller norm. Inductively, 
this gives the result. 

Suppose 7 is a loop of length n on T, with ^(7(2), c) > 2 and d(7(i±l), c) = 
(i(7(i),c) — 1. (That is, i is a position on the loop where 7 reaches a local 
maximum distance from the centre.) Consider the modified loop 7', which 
agrees with 7 except at position i, where it passes through the vertex 2 closer 
to the centre than 7(2) (possibly the central vertex itself). Consider also the 
'snipped' loop tc of length n — 2, obtained from 7 or 7' by removing the i-th 
and i + 1-th positions. We name the vertices as s = 7(2), r = 7(2 ± 1), and 
t = i{i). 

Applying a cap at position % to A, we have fl;(A) = 0. Evaluating this at 
n gives 

= y/dh^rf* a(A)(7r) = v/dim(s) fcl A(7) + y/dun{t) ki A^') . 

(Here ki is the number of critical points in the cap strand, either 1, 2 or 3 de- 
pending on the position of the point i around the boundary of the rectangular 
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box, as follows 



k = < 



1 when we have | ~| or ^ 

2 when we have ^ j j i or i 

3 when we have M or / 





although the case k = 3 never occurs for us as we always consider boxes with 
equal numbers of strands above and below.) 

For the second statement, consider some collapsed loop a of length n — 2. 
Write ctjj for the collapsed loop of length n which makes an extra visit to the 
j-th spoke of the graph at position i, and write Uj for the vertex adjacent to 
the central vertex on the j-th spoke. (If the graph has even supertransitivity, 
then i ought to be odd, while if the graph has odd supertransitivity, % is even.) 
Now, 



Using this formula, we can express the value of A on any collapsed loop which 
visits some spoke in terms of other collapsed loops which visit that spoke 
strictly fewer times. 

For the final statement, note that A has value zero on any loop which visits 
a univalent vertex at positions i and i + 2 by a similar argument as above. 
The value of A on a loop with 2k + 1 consecutive vertices either in a fixed 
spoke of length k — 1 or at the central vertex is then a multiple of the value of 
A on a loop that visits the end of that spoke twice consecutively by our first 
argument, and is thus also zero. □ 

Corollary A.2. If A is a lowest weight vector in a spherical graph planar 
algebra, and 7 denotes the 'collapsed' loop corresponding to 7, then 



The corresponding formula for a lowest weight vector in a lopsided graph 
planar algebra is 




(A.l) 




(A.2) 




(A.3) 
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where 



when we have 

1 when we have E 

















or 


1 1 






J 







2 when we have f 1 \. 



(The exponent is just the number of minima on the strand.) 

Rotation acts on the set of collapsed loops, so if we are trying to specify a 
lowest weight vector A which is also a rotational eigenvector, then it suffices 
to specify A only on a representative of each such orbit. 

Lemma A. 3. Fix u> an n-th root of unity. Suppose we have specified the 
values of an element A on a rotation representative of each collapsed loop 
which avoids the initial arm, and further that 

(a) if a representative is fixed by the k-fold rotation, and u k ^ 1, the corre- 
sponding value of A is zero, and 

(b) condition (3) of Lemma A.l holds where appropriate on these values, i.e., 
if a loop visits an arm of length k — 1 at least k consecutive times, the 
corresponding value of A is zero. 

Then we can make three consecutive well-defined extensions, defining, in turn, 
the values of A 

(1) on every collapsed loop avoiding the initial arm, using the condition that 
A is a rotational eigenvector with some eigenvalue u, 

(2) on every collapsed loop, via Equation (A.l), and finally 

(3) on every loop, via Equation (A. 2). 

The resulting element A is a rotational eigenvector, and is a lowest weight 
vector if and only if 



A(00...0) = (-!)» E ff[(^V l/2 M(^ 2 ---n) = 

Oi=l,...,m \i=l 



(here we denote the value of A on the collapsed loop which successively visits 
arms ai, a 2) . . . , a n by A(a\a2 ■ ■ ■ a n ) ). 
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Example A A. Consider the Haagerup principal graph, with n = 4 and three 
arms of length 3. There are 5 representatives under rotation of collapsed 
loops avoiding the initial arm, namely 1111,1112,1122,1212,1222 and 2222. 
Condition (3) of Lemma A.l says that A(llll) and A(2222) must be zero 
if we want to extend A to a lowest weight vector. Further, if uj — ±i then 
A(1212) must be zero also. Finally, the condition A(0000) = is trivially true 
when uj ^ 1, since there A(0000) = uA(0000), while it is non-trivial when 
to = 1. Thus the the rotational eigenspaces of lowest weight 4-boxes in the 
graph planar algebra are 3 dimensional for u = —i, and 4 dimensional for 
u = -1. 

Thus, in each of the following subsections, we list representatives of the 
rotational orbits of the collapsed loops avoiding the initial spoke. We then 
list the values of our generators on these loops, and refer to the function 
lowestWeightCondition in the Mathematica notebook Generators .nb (in- 
cluded with the arXiv sources of this article) for the elementary check of 
Lemma A. 3 that these actually determine a lowest weight vector with the de- 
sired rotational eigenvector. Note that this notebook takes quite a while to 
run on 4442, as it needs to lift the specified values, expressed as particular 
roots of their minimal polynomials, back to the fixed number field Q(/^4442) 
described below. This notebook also regenerates all the values in accordance 
with the above lemma, in a format compatible with the FusionAtlas package. 
This notebook, however, is completely independent of that package. We note 
that this method of describing lowest weight vectors was implicitly used in 
[BMPS09], but without explanation of why it is always possible. 



A. l 4442 

We give here two generators A and B which are rotational eigenvectors and 
lowest weight vectors, but are not self adjoint. We'll correct them in Subsection 

B. l by a phase to obtain self-adjoint elements. 

The two generators Aq, Bq for 4442 have rotational eigenvalues 

/ 3\ (2 

loa = exp 2ixi- cub = exp 2ixi 



5/ " V 5 

and we chose square roots a a, vb- 

( 8 \ / 2 \ 

a a = exp 2ni— as = exp 2ixi— . 



io J V 10 J 

We express A , B here by giving their coefficients on representatives of 
the rotational orbits of the collapsed loops. We write these coefficients as 
algebraic numbers, that is, roots of certain integer coefficient polynomials. 
Recall that the notation Ai„,...,ao indicates the root of the polynomial ^ a^x 1 
which is approximately equal to z. (We give z to sufficiently high precision 
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that it clearly distinguishes amongst the roots.) In fact, we know that all of 
these numbers lie in a single number field, but it is a terrifying one: Q(/X4442), 
where /Z4442 is the root of 

x 16 - 137624a; 15 + 8933996874a; 14 - 350479594607884a; 13 

+ 9011981487580477099a; 12 - 153965505437561352450336a; 11 

+ 1677614319697333636358399288a; 10 

- 10020910478354387137539071365292a; 9 

+ 14924720493503061057816631839615921a: 8 

+ 106599644350188183735570113321511839620a; 7 

+ 190438222993363250317220258290893662331668a; 6 

+ 186721586395717486269136749424837851231468948a; 5 

+ 116885475358299597752661285079387479235051171199a; 4 

+ 48661380081758389460116404892672194953761868153880a; 3 

+ 13292626196670539867585523753990034583320535147867050a; 2 

+ 2196889811116840388811767217587569058102181055950438988a; 

+ 171603625204099732635330174190840922237651443646904957481 

which is approximately 17589.4 + 13246. 7i. The polynomials in /Z4442 required 
to express these numbers themselves tend to be horrific (coefficients whose 
numerator and denominator may have hundreds of digits), and we chose to 
spare the reader from the danger of trying to read them. The overenthusi- 
astic may of course view them in the Mathematica notebook. Nevertheless, 
it is important to remember that the calculation of moments in §B must be 
performed inside this number field in order to be tractable. 

The reason we work with the non-self- adjoint generators Aq and Bq is 
simply that the coefficients of the phase corrected generators require an even 
larger number field; sufficiently large, in fact, that our computers can't perform 
the necessary calculations there! 



4)(mi2) 


= A; 


A (11113) 


= A 


^o(H122) 


= A; 


A (11123) 


= A, 


A (11132) 


= A 


A (11133) 


= A, 


A (11212) 


= A 


^o(H213) 


= A 


A (11222) 


= A; 


A)(H223) 


= A 


A (11232) 


= A 



(0.009547+0. 023000i) 

3941830656,370676736,12897792,-5087232,111600,36792,852,-6,1 
(0.026510+0.063867i) 

1247219856,379505736,55366092,5211378,385155,25407,1422,54,1 
(-0.37586+0.02986i) 

3941830656,2911168512,562837248,25208064,6989040,233784,7428,42,1 
(0.13774+0.1319H) 

77951241,-40286727,21922488,-4650534,645165,-27324,2988,-27,31 
(0.11518-0.15201i) 

1247219856,652215888,202145868,6602796,6554925,997056,165888,9768,496 
(0.22307-0.01772i) 

77951241,-34088769,8046702,-1582902,130005,-10098,2412,24,16 
(-0.15040+0.62850i) 

15397776,-13490712,3167748,-6096654,2141865,640836,414528,2883,961 
(-0. 57158-1. 17710i) 

1247219856,3177168624,5503926672,4667265522,2501734815,231250113,72236142,4203666,122461 
(0.024549-0. 102591i) 

3941830656,2911168512,562837248,25208064,6989040,233784,7428,42,1 
(-0.181705-0.437753i) 

1247219856,-598818096,213923592,-131814378,67404555,-2266407,34767,-279,1 
(-0. 23036+0. 96265i) 

1247219856,1315874160,1776121020,856453500,225182025,32508000,1964250,13125,625 
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A niTO^ — \(0-48303-0.03838i) 

Ji0(LlZ66) — -"\7795124i, -49106898,33170958,-32103216,9349425, -303966, 101448, 1392, 16 



A / lr !19\ - »(1.04250-0.79086i) 
A.o(L16lZ) — A 1 2 4 72 1 g 8 56 ! -884877696,195611112,-129102498,2093568525,564703758,25239672,-2431059,122461 

A (-\-\oio\ _ \ (-0.08426+0.35214i) 

^0l 1 - L <J i 'Jj — A 77951241,74852262,35328798, 11706039,2207655, 166914,12168,432,16 

(0.360178-0.308085i) 
V 35316, -24678,7389, 16 

(-0. 18246+0. 76250i) 



t4o(11322) — A353 1 (i. 2 Ki78.73S!).l(iS.] 



A0{IL6Z6) — ^77951241, -155425716,96700392, -79082163,68641425, -4964112,278712,-144,16 
A (-\-\ooo\ — \ (-0.41334-0. 25287i) 

Aoyi-L6d/,) — ^77951241, 121336947,113029263,60444549,17861715,2067174,73368,-1248,16 
A n 91 99\ - \(0.77203-0.06134i) 

Aq^LZLZZ) — ^15397776,-13490712,3167748,-6096654,2141865,640836,414528,2883,961 

A (I oioo\ _ \ (-0.3333+1.0259i) 
J\Q\i.ZLZH) — A 81, 108, 144,72, 16 

A M 91^9^1 - \ (-0.4911-0.9604i) 

^Ol+^ioZj — ^77951241, -75805794,99558072, -91700667,145832805, -61946748,27924912,-9145536,3041536 
4 M91^ — \(l-15226-0.09155i) 

^0(J-^4o^; — A7795124I, -151373205,49269465, 15230025, 11397915,918000,319500,6000,400 
A l-\O0-\-X\ — \ (-0-06347-0. 15291i) 

Aq{LZZ16) — A i247219856,-2526859800, 1974773520,-491005800,55811025, -10587375,2076750,-45000,625 

a /i 9999^1 — \ (0-075927+0. 046450i) 
A0{LZZZZ) — A 3941830656,370676736, 12897792, -5087232,111600,36792,852, -6,1 

A M999<^ - \(0.13497-0.38939i) 

^0{LZZZO) — A 7795i241,31228173, 13579083, 7130106, 2059425,324216,29448, 1458,31 
A n 99^ - \(0.49113-0.87531i) 

^Ol-L^^-J^J — A 1247219856,1670588064,-268120368,1191806622, 1842147765, -527829912,98070912, -6203229, 122461 
A M99^ _ x (-0.70611+0.05610i) 

^Oli^^-J-Jj — A 77951241,83195667,22443723, 14206509,7526925, -899316, 106308,-48,16 
A H9<I1^ — \(0-07491+1.38254i) 

JiQ{LZ6L6) — A 7795i241, -175449888,329852088,-416257056,333376560, -158331456,44336448, -6899328,476416 
A — \ (-0-55570+0. 83581i) 

Aq{IZ6ZZ) — A i247219856,-884877696,195611112, -129102498,2093568525,564703758,25239672, -2431059,122461 
A /19W\ — \ (-0.24343-0.67826i) 

A0{LZ6Z3) — A 7795i241,325154412,587305728,611166384,426416400,193325184,52769088, 7788672,476416 
A — \ (0-031564+0. 019310i) 

JiQ(LZ66Z) — A 7795i241, -134686395,185147775,-112077675,34665165,-5968350,586800,-24000,400 
A (-IQ-,o \ _ x (1.1727-0.5217*) 

J±0{ld±dZ) — ^190096,-298224,48892,140748,-31635,79392,9592,14424,29776 
A f-\Q-iQQ\ — \ (-l-3384+0.1063i) 

^o(lolOOJ — A ii881,31174,44152,52112,25680,11968,17152, -1024,4096 
A (TJ000\ - \(-0-26442+0.27576i) 

/i {LOZZZ) — ^3041536,-390656,219552,-53008,49325,-14677,3417,-389,31 
A f1399^ — \ (-0-90634-0. 55448i) 

Aq{16ZZ6) — A 190096,419432,437068,227614,74265, -12256, 2458,97,1 
A (-\1<yX0\ — \ (-0-6626+0. 0854i) 

JiQ{L6Z6Z) — A ig0096,-298224,48892,140748,-31635,79392,9592,14424,29776 
A (^1'y^1\ — \ (0-6692-0. 0532i) 

Aq(16Z66) — A n881, -15587,11038,-6514,1605, -374,268,8,16 
A — \ (-0-25174-0.60649i) 

Ji0[io6ZZ) — A i90096, 155216,93912,6998, -10335, -2908, 1022,39,1 
A — \(0-2574+0.6200i) 

Ji.0{L66Z6) — A H881, -3052, 1228, 1861,255, -1174,628, -152, 16 
a f 0000 o\ _ \ (0.026510+0. 063867i) 

7i0 [ZZZZ6) — A 1247 2i9856,379505736,55366092,5211378,385155,25407,1422,54,l 
A (090W\ - \ (0-22307-0.01772*) 

S±0{ZZZdd) — A 7795i241, -34088769,8046702, -1582902, 130005, -10098,2412,24,16 
A (99<iO-i\ - \(-0.08426+0.35214i) 

n.Q{ZZOZO) — ^77951241, 74852262,35328798, 11706039,2207655, 166914,12168,432,16 
A (Wl^ - \ (-l-3384+0.1063i) 

^0{ZOZOO) — A H881,31174,44152, 52112, 25680,11968, 17152, -1024,4096 
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B 
B 
B 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
Bo 
B 



(11112) 
(11113) 
(11122) 
(11123) 
(11132) 
(11133) 
(11212) 
(11213) 
(11222 
(11223 
(11232) 
(11233) 
(11312) 
(11313) 
(11322) 
(11323) 
(11332) 
(12122) 
(12123) 
(12132) 
(12133) 
(12213) 
(12222) 
(12223) 
(12232) 
(12233) 
(12313) 
(12322) 
(12323) 
(12332) 
(13132) 
(13133) 



x (0.075927+0.046450*) 

A 3941830656,370676736, 12897792 -5087232,111600,36792,852,-6,1 
x (-0.058988-0. 036087*) 

A 1247219856,379505736,55366092,5211378,385155,25407,1422,54,l 
x (0.024549-0.102591*) 

A 3941830656,2911168512,562837248,25208064,6989040,233784,7428,42,l 
x (0.11968+0.39436*) 

A 77951241, -40286727,21922488, -4650534,645165, -27324,2988, -27,31 
x (-0.28523+0.29746*) 

A 1247219856,652215888,202145868,6602796,6554925,997056, 165888,9768,496 
x (0.05208-0.21763*) 

A 77951241, -34088769,8046702 -1582902,130005,-10098,2412,24,16 
x (0.77203-0.06134*) 

A 15397776,-13490712,3167748 -6096654,2141865,640836,414528,2883,961 



A 

) = A 
) = A 

A 
A 
A 



(-0.68071-0.73758*) 

1247219856,3177168624,5503926672,4667265522,2501734815,231250113,72236142,4203666,122461 
(-0.37586+0.02986*) 

3941830656,2911168512,562837248,25208064,6989040,233784,7428,42,1 
(0.404307+0.247347*) 

1247219856,-598818096,213923592,-131814378,67404555,-2266407,34767,-279,1 
(-0.16504+0.01311*) 

1247219856,1315874160,1776121020,856453500,225182025,32508000,1964250,13125,625 
(-0.16485+0.68888*) 

77951241,-49106898,33170958,-32103216,9349425,-303966,101448,1392,16 
(-0.55570+0.83581*) 

1247219856,-884877696,195611112,-129102498,2093568525,564703758,25239672,-2431059,122461 
(-0.36094+0.02868*) 

77951241,74852262,35328798,11706039,2207655,166914,12168,432,16 

(0.360178-0.308085*) 
35316,-24678,7389,168,1 

(1.14251-0.09077*) 

77951241,-155425716,96700392,-79082163,68641425,-4964112,278712,-144,16 



x (-0.27155-0.65421*) 

A 77951241, 121336947,113029263,60444549,17861715, 2067174,73368, -1248,16 



(-0.15040+0.62850*) 

15397776,-13490712,3167748,-6096654,2141865,640836,414528,2883,961 



= A 
= A 



(-0.3333-1.0259*) 
81,108,144,72,16 

(0.7616-0.7639*) 

77951241,-75805794,99558072,-91700667,145832805,-61946748,27924912,-9145536,3041536 
(-0.008611+0.035986*) 

77951241,-151373205,49269465,15230025,11397915,918000,319500,6000,400 

(0.84435+0.51656*) 
1247219856, -2526859S 

(0.009547+0.023000*) 

3941830656,370676736,12897792,-5087232,111600,36792,852,-6,1 



x (0.84435+0.51656*) 

A 1247219856, -2526859800,1974773520,-491005800,55811025,-10587375,2076750,-45000,625 

A 



A 



(-0.18897-0.02576*) 

77951241,31228173,13579083,7130106,2059425,324216,29448,1458,31 



> (-1.29611+0.17986*) 

A 1247219856, 1670588064,-268120368, 1191806622, 1842147765, -527829912,98070912, -6203229, 122461 



(0.11277-0.47125*) 

77951241,83195667,22443723,14206509,7526925,-899316,106308,-48,16 
(0.56984-0.44112*) 



77951241,-175449888,329852088,-416257056,333376560,-158331456,44336448,-6899328,476416 

(1.04250-0.79086*) 
1247219856,-8848776 

(-1.29172+0.49847*) 

77951241,325154412,587305728,611166384,426416400,193325184,52769088,7788672,476416 

(0.44313+1.06757*) 
77951241,-1346863 
(-0.6626+0.0854*) 



> (1.04250-0.79086*) 

A 1247219856, -884877696, 195611112, -129102498,2093568525,564703758,25239672, -2431059,122461 

a; 

> (0.44313+1.06757*) 

A 77951241, -134686395, 185147775, -112077675,34665165 -5968350,586800,-24000,400 



190096,-298224,48892,140748,-31635,79392,9592,14424,29776 
(-0.3125+1.3058*) 

11881,31174,44152,52112,25680,11968,17152,-1024,4096 
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#0(13222) 


= a; 


B (13223) 


= a 


£ (13232) 


= a; 


S (13233) 


= a; 


B (13322) 


= A 


So (13323) 


= A 


S (22223) 


= A 


S (22233) 


= A 


S (22323) 


= A 


S (23233) 


= A 



(0.10678-0. 14092i) 

3041536,-390656,219552,-53008,49325,-14677,3417,-389,31 
(-0.27865-0.67130*) 

190096,419432,437068,227614,74265,-12256,2458,97,1 
(1. 1727-0. 5217i) 

190096,-298224,48892,140748,-31635,79392,9592,14424,29776 
(0.1562-0.6529i) 

11881,-15587,11038,-6514,1605,-374,268,8,16 
(-0.38319-0.23443i) 

190096,155216,93912,6998,-10335,-2908,1022,39,1 
(-0.5727-0.3503i) 

11881,-3052,1228,1861,255,-1174,628,-152,16 
(-0.058988-0.036087i) 

1247219856,379505736,55366092,5211378,385155,25407,1422,54,1 
(0.05208-0. 21763i) 

77951241,-34088769,8046702,-1582902,130005,-10098,2412,24,16 
(-0.36094+0.02868i) 

77951241,74852262,35328798,11706039,2207655,166914,12168,432,16 
(-0.3125+1.3058?) 

11881,31174,44152,52112,25680,11968,17152,-1024,4096 



A.2 3333 

The self-adjoint generators A, B of P 3333 have chiralities a A = 1 and 
Their values on collapsed loops are as follows: 



A(1112) = 


1 

4 


(2-V5) 


A(1113) = 


1(3^5-7) 


A(1122) = 


1 

4 


(V5-3) 


A(1123) = 




A(1132) = 


1 

8 


(3-V5) 


A(1133) = 




A(1212) = 


1 

4 


(VB-l) 


A(1213) = 




A(1222) = 


1 

4 


(2-V5) 


A(1223) = 




A(1232) = 


1 

8 


(VB-l) 


A(1233) = 


1 


A(1313) = 


1 

4 


(3-V5) 


A(1322) = 




A(1323) = 


1 

4 




A(1332) = 


1 


A(1333) = 


1 

8 


(3V5 - 7) 


A(2223) = 




A(2233) = 







A(2323) = 


1 

~2 


A(2333) = 


1 

8 


(3-V5) 






5(1112) = 


1 

8 


15) 


5(1113) = 




5(1122) = 


1 

2 


(2-V5) 


5(1123) = 


1 

~ 4 
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5(1132) = 


1 

4 


5(1133) = 


1 
2 


(2-\/5) 


5(1212) = 


I(5-*/S) 


5(1213) = 


1 

4 


(2-V5) 


5(1222) = 




5(1223) = 


1 

8 


^5 - 3^5) 


5(1232) = 




5(1233) = 


1 

8 


^3^5 - 7) 


5(1313) = 


1 

2 


5(1322) = 


1 

8 


^5 - 3^5) 


5(1323) = 


i(^- 5 ) 


5(1332) = 


1 

8 


^3\/5 - 7) 


5(1333) = 




5(2223) = 


1 

4 


^9 - 4>/5) 


5(2233) = 




5(2323) = 


1 

4 


J^3\/5 - 7) 


5(2333) = 


\ (9 - 4v£) 









Clearly all of these entries lie in the field Q(y/5). 



A. 3 3311 

The self-adjoint generators A, B of P 3311 have chiralities a a = i, 00 a = a\ = 



— 1, and (7b = cjb = 


= 1. Their values 


on collapsed loops are as follows 


A(1112) 


-i(»-5^) 


A(1113) = ^ (3V3-5) 


A(1123) 


= 1(^-0 


A(1132) = 1 (l - Vt) 


A(1212) 


= 1(5^-9) 


A(1213) = 1^-3) 


A(1232) 




A(1313) = - (7-3V3) 


A(1323) 




A(2323) = 



5(1112) = - (9-5V3) 
5(1123) = - (3- V^j 
5(1212) = ^ (l3V3-2l) 
5(1232) = - ^7>/3- 15) 
5(1323) = 1 (V3-l) 

Clearly these entries all lie in Q(\/3). 



5(1113) = -(^-3) 
5(1132) = - (3- V^j 
5(1213) = J 
5(1313) = ^ (l - V3) 
5(2323) = ^ 
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A.4 2221 



As in Subsection A.l, we work with non-self- adjoint generators for 2221. We 
will correct them by phases in Subsection B.4 to get self-adjoint elements. 
The generators Aq, Bo have rotational eigenvalues 

uja = exp ( 2ni- J ujb — 1 



for which we choose square roots 

a A = exp ^ 2? ™0 a B = 1. 

Their values on collapsed loops are as follows: 

A f-\-\o\ — \ (0-5762-0. 5412i) ■ f-i-io\_ \ (-0.4208+0.3953i) 

A.q{1LZ) — A 9)9j9) _ lgj21j _ 12)12) _ 61 A.q{116) — ^81,81,54,-27,-27,-9,6,3,1 

A M99\_ \ (0-3194+0. 0964i) . ,.„„>_ , (-0.6319-1. 0945i) 
JiQ{LZZ) — ^9,9,9,-18,21,-12,12,-6,1 — A 9,9,12,-3,l 

A C\o<y\ — \ (-l-1319+0.5621i) , fOOQN _ x (-0-4208+0. 3953i) 

JiQ{l6Z) — A 9)27)30j12j1 Aq{ZZ6) — ^81,81,54,-27,-27,-9,6,3,1 



S (112) 
B„(122) 
Bo(132) 



(-0.0605-0.1502i) 
225,-45,288,-45,201,-39,36,3,1 

(0.3396+0.8433i) 
225,-45,288,-45,201,-39,36,3,1 

3 + \/21 N 



These entries all lie in the number field 

a; 8 - 18x 7 - 345x 6 



So(H3) 
B (123) 

So (223) 



A 



(-0.3528-0.8759i) 
45,36,42,3,-1 

(-0.9110+0.8759*) 
45,54,24,-39,5 

(-0.3528-0.8759i) 



45,36,42,3,-1 

2221 ) where /Z2221 is the root of 



7146x 5 + 84726x 4 - 1458918x 3 



13821786x 2 + 101759328a; + 1245393549 



which is approximately 14.85 + 9.90i. 

In fact, we could have tried working directly with self-adjoint generators, 
but this would have required a degree 32 number field. We could still calculate 
the requisite moments; however, we could not express the values of the gen- 
erators on collapsed loops in the above compact form, and instead, we would 
have to write out the 32 coefficients in the number field for each value. 



B Moments 

In the following subsections, we give the quadratic and cubic moments for our 
generators calculated directly from the graph planar algebra. Recall that the 
generators given in Appendix A are not normalized, so Tr(A 2 ) and Tr(B 2 ) 
need not be equal to 1. 
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B.l 4442 



Recall from Appendix A.l that the generators given for 4442 are not self- 
adjoint, since it would require working in an even larger number field (as if 
Q(/i4442) isn't terrifying enough)! However, for each non self-adjoint So £ %$o, 
there is a 7s E U(l) such that S = -fsSo is self-adjoint. Hence we compute the 
moments in the graph planar algebra using the non self-adjoint generators, 
and we correct them afterward to agree with what we would get from first 
making our generators self-adjoint. This amounts to multiplying the moment 
by the correction factor for each generator that appears, e.g., 

Tt(A 2 B) = Ti{{ lA A ) 2 lB B ) = "fan T±(A*B ) . 

In the case of 4442, the correction factors for Ao,Bq are equal and given by 



lA = IB 



'437 135^ 1 . /15 



872 872 
and the corrected moments are as follows: 



+ 436 2 Vt( 15943 + 3933v/K )' 



Tr(v4 2 ) = 
Tr(AB) = 



12706131015 5682354525\/5 
11881 11881 



2> ( , 12706131015 5682354525\/5 
' 1 J ' ,| ; 11881 + 11881 



. :! ,'36768021311025 16443159009075\/5 
' ' ' ' ^ 2590058 + 2590058 

/25201176715665 11270308849830^ 

Tt(A B) = -\ + — 

1 1295029 1295029 

/25201176715665 11270308849830^ 

Tt(AB ) = -\ + — 

1 1295029 1295029 

p:u / 36768021311025 16443159009075^ 
' { J ' ,|; 2590058 + 2590058 



,12706131015 5682354525\/5 
Tr (^ 2 = \ 110 ci + 



Tr(AB) = 
Tr^ 2 ) = 



11881 11881 

12706131015 5682354525^5 
11881 11881 
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, : . ,'36768021311025 16443159009075^ 
1 ' ' A ' " l " 2590058 + 2590058 

/25201176715665 11270308849830^/5 

Tr A Z B) = -\ + — 

1295029 1295029 

/ 25201176715665 11270308849830\/5 
' ' ' 4 J ' ~ V 1295029 + 1295029 

,'36768021311025 16443159009075^ 
Tr (B*) = \ — — + 



2590058 2590058 



B.2 3333 



Since the generators of 3333 are self-adjoint, no corrections are needed. The 
moments are as follows: 

Tr(A 2 ) = X - (ll + 5v^) Tr(i 2 ) = ^ (ll + 5v^) 

Tr(AB) = Tr(AB) = 

Tr(B 2 ) = 6 + 3^5 Tr(B 2 ) = 6 + 3^5 

Tr(A 3 ) = ^ (-27 - 12-v^) Tr(i 3 ) = 



Tr(A 2 5) = 1 (-25 - 11 Vs) Tr(A 2 5) = -^-f + ^ 

Tr(AB 2 ) = (63 + 27v^) Tr(AB 2 ) = 

Tr(B 3 ) = i (l5 + 6v^) Tr(fi 3 ) = 3^ (9 + 4v^) 
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B.3 3311 



Since the generators of 3311 are self-adjoint, no corrections are needed. The 
moments are as follows: 



Tr(A 2 



12 + 8^3 



Tr(i 2 



1 



12 + 8^3 



Tr(AB) = 



Tr(5 2 ) 
Ty(A 3 ) 

Ti(A 2 B) 

Ty{AB 2 ) 

Tr(5 3 ) 

B.4 2221 



1 



3 
] 

27 

1 

9 



9 + 4^) 

I (-3 - 2V3 
1 (15 + 14V3 
i (23 - 2V3 



Tr(AB) = 
Tr( J B 2 ) = ^(9 + 4v / 3) 
Tr(i 3 ) = 



Ti(A 2 B) 
Tt{AB 2 ) 
Tt(B 3 ) 



'44 76 
+ 



27 27V3 



'9 



4 4^3 



As in Subsection B.l for 4442, we now correct our non-self-adjoint generators 
for 2221 by multiplying by the following phases: 



1A 



IB 



1 -V21 



19 V21 1 
+ - i\ - 9 + 19V21 

20 20 10 V 2 



which yield the corrected moments: 



Tr(A 2 ) = 3\/55 I2\/2I 
Tr(AB) = 



Tr(A 2 ) = 3\/55 |2\/2I 
Tr(AB) = 



Tr(5 2 ) 
Tr(A 3 ) 

Tr(A 2 5) 
Tt(AB 2 ) 

Tr(B 3 ) 



'60093 13113V21 



50 

-18-4\/21 



50 



'4413 963 V21 
+ 



10 



10 



198 + 43V21 



11(5(57 25 17 v 21 



250 



250 



Tr(f? 2 ) 
Tr(i 3 ) 

Tr(A 2 B) 
Tr (AB 2 ) 

Tr(£? 3 ) 



(50093 . 131 1 3v/2l 



50 

-18 -Ay/21 



50 



4413 963V21 
+ 



10 



10 



198 + 43V21 



11667 2547V21 
+ 



250 



250 
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